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Abstract 

Successful  simulation  of  kinetics  of  small  fatigue  crack  growth  entails  three  aspects:  Stage  I,  Stage  II  growth  rate 
prediction  and  transition  prediction.  In  this  paper  attention  is  focused  on  growth  rate  predictions.  By  using 
microstructurally-affected-zone  and  process  zone  concepts,  microscopic  fatigue  behaviour  of  small  fatigue  crack 
propagation  is  logically  linked  with  macroscopic  fatigue  behaviour,  showing  an  intrinsic  relation  between  small 
fatigue  crack  growth  and  macroscopic  low-cycle  fatigue  properties  during  crack  growth.  Furthermore,  variation  of 
relatively  big  plastic  zone  size  associated  with  a  growing  small  fatigue  crack  is  kinetically  simulated.  As  a  result  a 
quantitative  prediction  model  of  growth  rates  for  Stage  I  and  Stage  II  growth  has  been  developed  whose  explicit 
advantage  is  that  the  growth  rate  of  small  fatigue  crack  can  now  be  predicted  in  terms  of  bulk  fatigue  properties  in 
conjunction  with  local  microstructural  characteristics. 


1.  Introduction 

Small  fatigue  crack  in  aluminium  alloys  has 
been  well  classified  into  microstructurally  small 
fatigue  crack  and  physically  small  fatigue  crack. 
Small  fatigue  cracks  may  initiate  directly  either  at 
slip  bands  in  matrix  as  a  result  of  dislocation 
motion  or  at  a  particle  bonded  in  matrix  due  to 
stress  concentration.  Wherever  a  small  fatigue 
crack  initiates  it  may  experience  two  stages  of 
growth  [1].  In  Stage  I  growth  there  are  two 
prominent  features  associated  with  the  small  fa¬ 
tigue  crack  growth:  shear  mechanism  primarily 
dominates  the  propagation  of  the  small  crack; 
and  local  microstructure  has  a  significant  effect 
on  the  growth  of  small  fatigue  crack.  This  influ- 
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ence  of  local  microstructure  cannot  be  negligible 
until  the  growth  mode  of  shear  mechanism  com¬ 
pletely  shifts  to  tension  mechanism.  In  Stage  II 
growth  small  fatigue  crack  growth  is  still  affected 
to  some  degree  by  the  local  microstructure.  How¬ 
ever,  the  propagation  of  small  fatigue  crack  is 
considered  to  be  different;  it  is  primarily  gov¬ 
erned  by  tension  mechanism  and  there  is  no 
longer  a  strong  dependence  of  propagation  upon 
the  local  microstructure  because  the  metallurgi¬ 
cal  influence  of  local  microstructure  on  fatigue 
crack  growth  is  averaged.  Thus,  the  mechanics  of 
continuum  may  be  applicable.  Furthermore,  with 
the  increase  of  small  fatigue  crack  length  up  to 
such  a  degree  free  from  influence  of  surface 
effect,  the  growth  is  inevitable  to  be  in  compli¬ 
ance  with  the  bulk  fatigue  behaviour.  This  implies 
that  there  exists  a  coherent  connection  of  fatigue 
growth  with  bulk  fatigue  properties,  in  other 
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words,  a  connection  between  micro-growth  of 
small  fatigue  crack  and  macro-fatigue  behaviour. 

Although  there  have  been  a  number  of  models 
proposed  to  predict  growth  rate  of  small  fatigue 
crack,  their  general  application  to  fatigue  life 
prediction  requires  the  determination  of  an  ex¬ 
pression  for  crack  growth,  which  is  obtained  from 
careful  and  tedious  characterisation  of  small  fa¬ 
tigue  cracks.  Particularly,  the  metallurgical  char¬ 
acteristics  of  polycrystalline  material  in  relation 
to  small  fatigue  crack  growth  has  not  been  mod¬ 
elled  yet.  The  coherent  connection  of  micro¬ 
growth  of  small  fatigue  crack  with  macro-fatigue 
behaviour  is  still  something  unclear.  As  a  result 
of  theoretical  study  and  experimental  observa¬ 
tions,  we  are  trying  to  physically  model  a  kinetic 
relation  between  metallurgical  characteristics  of 
material  and  characteristics  of  small  fatigue  crack 
tip  during  crack  growth,  and  hence  to  make  a 
connection  of  micro-fatigue  with  macro-fatigue. 
A  kinetic  variation  of  plastic  zone  in  the  process 
of  fatigue  crack  propagation  is  also  simulated 
with  emphasis  on  the  variation  of  the  ratio  of 
plastic  zone  size  to  small  crack  length.  This  simu¬ 
lation  is  totally  based  on  experimental  results  of 
plastic  zone  measurement  along  with  elastic-plas¬ 
tic  fracture  mechanics.  In  general,  we  are  trying 
to  establish  a  general  model  for  aluminium  alloys 
that  can  predict  growth  rate  by  means  of  conven¬ 
tional  low-cycle  fatigue  properties  without  te¬ 
dious  crack  length  measurements  and  can  simu¬ 
late  the  kinetic  evolution  of  small  fatigue  crack 
growth. 


2.  Background  of  modelling 

2.1.  Simulation  of  microstructure  for  polycrystalline 
materials 

Microstructure  has  a  dramatic  influence  on 
small  fatigue  crack,  particularly  on  its  propaga¬ 
tion  of  early  stage.  Since  a  growing  small  fatigue 
crack  is  controlled  to  some  extent  by  local  mi¬ 
crostructure,  it  seems  to  us  that  local  microstruc¬ 
ture  should  be  modelled.  However,  it  is  impossi¬ 
ble  to  simulate  each  local  effect  of  single  metal¬ 
lurgical  feature,  such  as  grain  size,  grain  orienta¬ 
tion  etc.,  on  propagation  of  small  fatigue  crack. 
What  is  known  is  that  it  is  possible  to  take  the 
microstructure  effect  into  account  by  defining  an 
appropriate  microstructural  dimension.  Proposed 
is  a  concept  of  microstructurally-affected-zone 
that  is  in  front  of  a  small  fatigue  crack  tip  to 
represent  local  metallurgical  features.  Its  size  can 
be  taken  as  an  appropriate  microstructural  di¬ 
mension  to  measure  the  effect  of  local  mi¬ 
crostructure.  Now  the  question  is  how  to  relate 
this  microstructurally-affected-zone  to  a  physical 
quantity  and  determine  it  quantitatively. 

Based  on  a  microscale  observation  of  a  real 
material,  it  was  stated  in  [1]  that  “the  microstruc¬ 
tural  features  in  metals  cause  break  up  of  the 
crack  front  into  segments  that  are  related  to 
elementary  blocks  operating  with  some  degree  of 
independence  from  their  neighbours  but  under 
the  general  influence  of  macroscopic  crack  of 
which  they  are  a  part”.  Following  this  concept  a 


Fig.  1.  Schematic  representation  of  (a)  elementary  block  and  (b)  tip  characteristic  of  a  small  fatigue  crack. 
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real  aluminium  alloy  is  conceived  to  be  composed 
of  elementary  blocks  of  finite  linear  dimension 
p* .  The  elementary  blocks  divide  the  material 
into  segments  and  are  identical  in  size.  Since  a 
surface  small  fatigue  crack  propagates  across  ele¬ 
mentary  blocks,  we  can  naturally  and  reasonably 
associate  the  elementary  block  size  with  the  mi- 
crostructurally-affected-zone  size  as  if  a  real  alu¬ 
minium  alloy  would  be  zoned  by  microstruc- 
turally-affected-zone.  Furthermore,  the  mi- 
crostructurally-affected-zone  size  is  taken  as  a 
measure  to  be  in  accordance  with  the  elementary 
block  size.  Therefore,  any  effect  of  local  mi¬ 
crostructure  on  a  growing  small  fatigue  crack  is 
thought  to  be  reachable  via  each  elementary 
block/microstructurally-affected-zone.  With  such 
a  simulation  we  are  able  to  equate  the  mi- 
crostructurally-affected-zone  size  to  a  simplified 
physical  quantity,  i.e.,  the  elementary  block  size. 
Fig.  1  provides  a  schematic  representation  of 
elementary  blocks  along  with  crack  tip  configura¬ 
tion. 

2.2.  Representation  of  small  fatigue  crack  tip 

Consider  a  physically  small  fatigue  crack  that 
lies  in  an  infinite  body  subjected  to  a  uniform 
tension  stress  (see  Appendix  A).  No  matter  how 
small  a  crack  tip  is,  it  is  still  valid  to  characterise 
a  small  crack  tip  as  surrounded  by  a  plastic  zone 
and  a  so  called  process  zone  just  ahead  of  the 
crack  tip  to  avoid  stress/strain  singularity  there. 
The  process  zone  is  physically  defined  to  be  de¬ 
pendent  upon  the  plastically  deformed  zone  and 
local  microstructure,  such  as  grain  size.  Naturally, 
the  process  zone  size  d  *  and  the  plastic  zone  size 
rp  are  mostly  concerned. 

However,  the  small  fatigue  crack  tip  compared 
with  relative  large  elementary  blocks  is  still 
smaller.  The  behaviour  of  crack  tip  is  believed  to 
be  controlled  by  localised  driving  force  in  a 
smaller  zone  at  the  crack  tip.  The  stress/strain 
distribution  varies  throughout  cracked  blocks  so 
the  uniform  stress/strain  contributing  to  lo¬ 
calised  driving  force  for  small  fatigue  crack  prop¬ 
agation  is  considered  to  occur  only  within  the 
process  zone  just  ahead  of  the  crack  tip.  This 
assumption  is  obviously  objective  because  the 


process  zone  is  defined  at  a  small  fatigue  crack 
tip  within  a  very  small  region  (smaller  than  the 
plastic  zone)  to  avoid  stress  singularity  at  the 
crack  tip.  The  process  zone  size  d  *  varies,  in¬ 
creasing  all  the  way  in  the  process  of  small  fa¬ 
tigue  crack  propagation  as  the  plastic  zone  size 
does.  An  incremental  advance  of  length  d*  is 
associated  with  a  local  measure  of  fatigue  dam¬ 
age  that  controls  fatigue  crack  growth.  Since  the 
total  increment  d  *  of  crack  length  is  a  direct 
result  of  fatigue  cycles  necessary  to  fracture  ma¬ 
terial  of  size  d  *  in  front  of  the  growing  small 
fatigue  crack,  the  crack  growth  rate  at  any  time 
will  be  therefore  given  by  the  process  zone  size 
divided  by  its  fatigue  damage  “life”,  which  can  be 
measured  quantitatively  by  the  ratio 
da  d * 


3.  Description  of  modelling  process  for  stage  II 
small  fatigue  crack 

It  has  been  modelled  already  that  small  fatigue 
crack  propagation  is  the  result  of  successive 
cracking  over  the  distance  d  * .  This  cracking  pro¬ 
cess  can  be  regarded  as  strain  controlled  fatigue 
in  the  process  zone  and  hence  the  fatigue  cycles 
can  be  determined  in  terms  of  Coffin-Manson 
relation.  In  this  respect  local  plastic  strain  in  the 
process  zone  significantly  contributes  to  the  driv¬ 
ing  force  for  crack  growth.  It  is  the  localised 
inelastic  strain  that  actually  controls  small  fatigue 
crack  propagation  and,  therefore  the  fatigue  life. 
In  addition,  most  polycrystalline  materials  can  be 
treated  as  isotropic  continuum.  Therefore,  free 
surface  effect  due  to  local  unusual  isotropic  (or 
anisotropic)  dislocations  [2]  on  early  stage  growth 
of  small  fatigue  crack  is  put  aside  in  the  mod¬ 
elling  for  the  time  being. 

3.1.  Strain  calculation  at  small  fatigue  crack  tip 

Consider  a  microcrack  of  length  2  a  that  lies  in 
an  infinite  body  subjected  to  uniform  tension 
loading.  The  microcrack  has  a  similar  configura- 
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tion  of  quasi-ellipse  whose  semi-major  axis  takes 
value  a.  Because  of  its  being  small,  the  microc¬ 
rack  can  be  mathematically  treated  as  a  “line 
crack”.  The  displacement  of  this  line  crack  in  an 
infinite  sheet  under  uniaxial  tension  stress  a  are 
given  for  plane  strain  condition  by  the  relation  [3] 


where  a  represents  half  surface  crack  length  and 
6  is  the  angle  between  the  point  of  interest  and 
the  major  axis  of  the  ellipse.  After  a  forward 
loading  the  microcrack  opens  becoming  a  mi¬ 
cronotch  with  blunted  tip  of  radius  d*.  In  the 
meanwhile,  a  plastic  zone  of  size  rp  is  formed  as 
a  result  of  applied  loading  at  the  micronotch  tip. 
The  opening  displacement,  8,  corresponding  to 
this  blunted  surface  micronotch  can  be  evaluated 
(See  Appendix  A)  by  the  following  expression: 


2(1  _  v2) 


(fl  +  rp)  1- 


(3) 


To  determine  local  strain  at  the  microcrack  tip 
we  adopt  a  physical  interpretation  of  crack  tip 
blunting  process  made  by  others  [4-6].  The  mi¬ 
crocrack  tip  blunts  to  a  semi-circle  of  radius  d  * 
without  singularity  of  local  strain  at  the  tip.  It  has 
been  theoretically  established  and  practically  em¬ 
ployed  [7,8]  that  the  normal  strain  distribution  in 
front  of  the  microcrack  tip  along  the  major-axis 
takes  a  form  of: 


(*>«), 


(4) 


where  17  is  the  deformation  in  a  yielded  zone  just 
at  the  crack  tip.  In  order  to  avoid  strain  singular¬ 
ity  at  x  =  a,  we  adopt  an  assumption  as  men¬ 
tioned  already  that  the  deformation  is  held  con¬ 
stant  within  the  distance  d  *  from  tip  of  the 
microcrack  that  has  been  named  “process  zone 
size”.  If  an  elementary  block  of  length  p ’  is 
initially  elongated  by  8/2  in  the  forward  loading, 
we  should  have 


(5) 


where  w  is  a  triaxiality  factor  taking  a  value  of 
2.57  [9].  To  meet  the  needs  of  continuity  of  Eq. 
(4)  with  relation  to  Eq.  (5)  at  x  =  a  +  d  * ,  77  must 
equal  8d  */(2p*«)  leading  to  the  following  strain 
expression: 


8d * 

2p*<o(x  —  a)  ’ 


(6) 


which  scales  strain  distribution  in  front  of  the 
microcrack.  Substituting  Eq.  (3)  into  Eq.  (6),  it 
yields  an  expression  of  strain  distribution  in  front 
of  the  elliptical  microcrack  tip  in  the  form  of 


d'(l-v2) 

p*<o£(x-a) 


X 

2' 

--1 

a 

1  +  — 

a 

At  x  =  a  +  rp,  it  is  a  boundary  between  local 
elastic-plastic  region  and  elastic  region  outside. 
The  strain  at  this  point  can  be  determined  by 
replacing  x  with  (a  +  rp)  in  Eq.  (7).  With  cyclic 
stress  imposed  on  the  infinite  body,  it  is  assumed 
that  the  cyclic  strain  distribution  can  be  evaluated 
by  means  of  Eq.  (7)  by  replacing  e,  cr  with  Ae/2 
and  Ao/2. 


3.2.  Kinetic  simulation  of  plastic  zone  associated 
with  a  growing  small  fatigue  crack 

Plastic  zone  size  associated  with  a  small  fa¬ 
tigue  crack  has  been  found  to  be  another  impor¬ 
tant  physical  quantity  that  plays  an  essential  con¬ 
trolling  role  in  small  fatigue  crack  growth.  Plastic 
deformation  at  a  small  crack  tip  dominates  the 
local  driving  force  around  the  crack  tip.  Unlike 
that  of  a  long  crack,  the  development  of  plastic 
zone  of  a  small  fatigue  crack  is  very  complicated 
due  to  local  microstructure.  It  has  been  found 
that  the  plastic  zone  is  likely  to  be  truncated  by 
grain  boundaries  due  primarily  to  misorientation 
between  neighbour  grains  when  the  small  fatigue 
crack  approaches  a  grain  boundary  [10,11].  A 
recent  experimental  measurement  [10]  on  7xxx 
aluminium  alloys  of  plastic  zone  indicated  that 
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the  relatively  large  plastic  zone  corresponds  to 
smaller  crack  length  and  the  ratio  of  plastic  zone 
size  to  crack  length  rp/a  initially  decrease  rapidly 
with  increase  of  crack  length  and  then  decreases 
gradually  until  approaching  a  relatively  stable 
value.  This  evolution  follows  a  power  function. 
On  the  other  hand,  it  has  been  recognised  that 
high  ratio  of  applied  stress  to  yield  stress,  <rmax/<ry, 
has  a  positive  contribution  to  the  comparatively 
large  value  of  the  ratio  rp/a  [12,13]  although 
some  deforming  grains  are  not  entirely  plastically 
deformed.  Based  on  previous  experimental  mea¬ 
surements  of  plastic  zone  size  as  shown  in  litera¬ 
ture  [10,11],  the  evolution  of  this  process  with 
respect  to  small  crack  length  and  stress  level  is 
illustrated  in  Fig.  2  which  shows  a  general  de¬ 
creasing  trend  of  ratio  rp/a  with  increase  of 
small  fatigue  crack  length.  Despite  the  fact  that 
the  plastic  zone  size  corresponding  to  a  given 
crack  length  may  change  in  different  aluminium 
alloys,  it  is  believed  that  the  evolution  of  the 
normalised  plastic  zone  size  rp/a  with  respect  to 
the  crack  length  and  the  stress  level  may  be 
approximately  identical.  However,  this  evolution 
of  normalised  plastic  zone  size  is  questionable 
when  it  is  extended  to  very  short  crack  size.  This 
point  is  going  to  be  addressed  later. 

As  has  been  pointed  out,  a  growing  small  crack 


Surface  Crack  Length,  2a  (mm) 

Fig.  2.  An  experimentally-determined  evolution  [10,11]  of 
normalised  plastic  zone  sizes  for  a  growing  small  fatigue 
crack. 


Surface  Crack  Length,  2a  (mm) 

Fig.  3.  Evolution  of  the  normalised  process  zone  size  in 
conjunction  with  small  fatigue  crack  length  and  microstruc- 


is  directly  associated  to  the  process  zone  just 
ahead  of  its  tip.  The  defined  process  zone  pro¬ 
vides  a  linkage  of  blunting  nature  of  the  crack  tip 
with  each  local  elementary  block  but  its  size  is 
under  a  general  influence  of  plastic  deformation 
at  the  microcrack  tip.  Once  a  small  crack  grows 
or  its  tip  deflects,  the  process  zone  size  changes 
simultaneously  in  compliance  with  local  mi¬ 
crostructure  constraint,  depending  upon  local 
plastic  deformation  and  local  microstructure.  Ow¬ 
ing  to  the  known  analytical  expressions  of  Eqs. 
(2)  to  (7),  use  can  be  made  of  the  requirement  for 
strain  equilibrium  at  the  boundary  between  the 
elastic-plastic  zone  ahead  of  the  microcrack  tip 
and  the  elastic  region  outside.  The  strain  at  the 
boundary  can  be  determined  in  terms  of  Eq.  (7). 
The  elastic  strain  just  outside  the  elastic-plastic 
zone  conforms  to  elastic  field  of  the  infinite  body 
which  can  be  directly  derived  in  terms  of  partial 
differentiating  with  respect  to  Eq.  (2),  i.e.,  Cdv/by) 
on  the  condition  that  the  microcrack  is  so  small 
that  it  can  be  regarded  as  a  line  crack.  Conse¬ 
quently,  we  derived  an  analytical  expression  that 
revealed  the  evolution  of  this  process  zone  size 
for  a  growing  small  crack  as  follows: 
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This  analytical  expression  shows  a  coherent  rela¬ 
tion  of  process  zone  size  d  *  with  microstruc- 
turally-affected-zone  size  p*  and  the  normalised 
plastic  zone  size  rp/a  via  stress  ratio  o-max/o-y 
that  makes  it  possible  to  link  the  microscopic 
fatigue  behaviour  of  small  fatigue  crack  propaga¬ 
tion  to  macroscopic  fatigue  behaviour  of  bulk. 
The  evolution  of  a  normalised  process  zone  size 
d*/p*  is  depicted  in  Fig.  3.  It  is  also  clear  that 
the  evolution  of  the  normalised  process  zone  size 
d*  /p  *  illustrated  in  Fig.  3  has  a  similar  pattern 
of  that  of  the  normalised  plastic  zone  size  rp/a 
shown  in  Fig.  2.  This  consistence  of  pattern  be¬ 
tween  evolution  of  the  normalised  plastic  zone 
size  rp/a  and  that  of  the  normalised  process  zone 
size  d*/p*  is  not  an  occasional  similarity  but 
reveals  and  provides  a  mechanical  basis  indicat¬ 
ing  that  the  process  zone  size  may  be  taken  as  a 
physical  quantity  to  kinetically  describe  small  fa¬ 
tigue  crack,  and  hence  is  correlated  with  a  small 
crack  growth  rate. 


3.3.  Coherent  connection  of  growth  rate  with 
macroscopic  fatigue  properties 


Successive  cracking  of  the  process  zone  over 
each  elementary  block  results  in  crack  extension. 
The  inelastic  strain  within  this  process  zone  takes 
a  primary  contribution  to  local  driving  force  at  a 
small  crack  tip.  Due  to  the  localised  nature  of 
crack  tip,  the  inelastic  strain  in  the  process  zone 
is  taken  to  be  larger,  comparing  with  elastic  strain. 
The  contribution  of  elastic  strain  can  be  ignored 
i.e.  Aep  ®  Ac.  Now  the  local  cyclic  inelastic  strain 
in  the  process  zone  can  be  theoretically  deter¬ 
mined  by  means  of  Eq.  (7),  letting  x  =  a  +  d* . 
This  yields  a  microstructurally-related  expression 
for  local  cyclic  inelastic  strain  at  the  small  fatigue 
crack  tip 


From  the  view  point  of  continuum,  any  influence 
of  local  metallurgical  features  on  small  fatigue 
crack  is  not  isolated  to  individually  take  effect  but 


generally  is  related  in  nature  to  macroscopic  fa¬ 
tigue  behaviour.  This  point  lays  a  foundation  in 
our  modelling  which  makes  it  possible  to  link 
physically  small  fatigue  crack  growth  to  bulk  fa¬ 
tigue  behaviour.  Since  the  failure  of  material 
element  at  crack  tip  is  essentially  due  to  strain- 
controlled  fatigue,  the  fatigue  cycles  necessary  to 
rupture  the  material  element  can  be  determined 
using  conventional  Coffin-Manson  relation  be¬ 
cause  relatively  large  plastic  deformation  does 
occur  at  a  small  fatigue  crack  tip  [10-14],  The 
cyclic  inelastic  stress-strain  response  of  alu¬ 
minium  alloys  obeys  a  well-familiar  power  func¬ 
tion,  A<r/2  =  /C'(Aep/2)n',  in  cyclic  fatigue.  A 
quantitative  relation  can  be  deduced  to  link  fa¬ 
tigue  cycles  to  failure  to  the  local  inelastic  cyclic 
strain  at  the  small  crack  tip  in  the  form  of 


Ae, 

2 


(  Of  -  °m)£'f 
l  K' 


1/1  +n' 

|  (2N{)b+c/1+n ,  (10) 


where  E  is  the  Young  modulus  (MPa);  K'  the 
cyclic  strength  coefficient  (MPa);  n’  the  cyclic 
strain  hardening  exponent;  Nf  the  fatigue  cycles 
to  failure;  b  and  c  the  material  constants  involved 
in  Coffin-Manson  relation  and  o-m,  o-{,  e'f  stand 
for  the  mean  stress  level,  fracture  stress  and 
fracture  strain  respectively.  Making  use  of  Eqs. 
(1),  (9),  (10),  growth  rate  of  a  physically  small 
fatigue  crack  can  be  given  as 

da 

d 77 


(11) 


This  is  a  growth  law  equation  for  a  single  physi¬ 
cally  small  fatigue  crack  by  which  growth  rate  is 
predicted  using  conventional  low-cycle  fatigue 
properties  and  two  proposed  microscopic  quanti¬ 
ties,  i.e.,  microstructuraily-affected-zone  size  and 
process  zone  size.  When  stress  level  in  Eq.  (9) 
approaches  fatigue  limit  (Acr)FL,  Eq.  (11)  pre¬ 
dicts  a  minimum  growth  rate.  This  equation  is 
applicable  to  small  fatigue  cracks  in  Stage  II 
propagation. 
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4.  Modelling  of  Stage  I  growth  rate 

Classified  in  [1]  are  two  stages  of  fatigue  cracks 
growth  no  matter  whether  the  crack  initiates  at 
slip  bands  or  emanates  from  a  particle.  It  was 
suggested  that  Stage  I  growth  is  similar  to  the 
process  by  which  a  crack  is  nucleated  in  an  active 
slip  band  through  the  formation  of  a  groove  by 
unequal  amounts  of  forward  and  reverse  slip  in 
neighbouring  packets  of  slip  planes.  In  Stage  I  a 
microcrack  initiates  and  propagates  in  a  crystallo¬ 
graphic  mode.  Shear  mode  growth  of  the  microc¬ 
rack  is  attributed  to  the  local  damage  due  to 
strongly  discrete  glide  edge  dislocations  at  the 
crack  tip.  Crack  propagation  is  governed  by  cyclic 
shear  stress  amplitude  At.  Barriers,  such  as  grain 
boundaries,  produce  perturbations  on  the  growth 
rate  of  small  crack.  Owing  to  a  complex  stress 
situation  around  a  small  fatigue  crack  tip,  Stage  I 
growth  features  a  mixed  mode  propagation.  Shear 
mechanism,  however,  is  a  primary  one  responsi¬ 
ble  for  the  Stage  I  growth.  A  theoretically  devel¬ 
oped  alternative  model  of  the  blocking  of  disloca¬ 
tions  at  grain  boundaries  [15,16]  provided  a  use¬ 
ful  solution  of  plastic  displacement  to  simulate 
the  discontinuous  character  of  slip  transfer  across 
grain  boundaries.  For  edge  dislocations  and  plane 
strain  condition,  the  shear  displacement  at  crack 
tip  is  given  [15,16]  by 

4(1  -  p2)  y/l-n2 

A  5  =  -— — 1 - Ara,  (12) 

E  n 

where  n  is  a  normalised  plastic  zone  size  defined 
by  n  =  a/(a  +  rp)  =  a/a0  for  Stage  I  propaga¬ 
tion.  Once  the  plastic  zone  extends  to  each  grain 
boundary  the  critical  value  nc  can  be  determined 
by  [16] 


where  Tmmp  is  a  comparison  stress,  rLi  =  tfl/V7 
O'  =  1,  3,  5,...)  and  tfl  is  a  shear  stress  corre¬ 
sponding  to  fatigue  limit  as  specified  in  [15,16]. 

Eq.  (12)  predicts  a  series  of  drops  of  growth 
rate  because  of  barrier  of  grain  boundaries  [16]  in 
the  process  of  small  fatigue  crack  growth.  How¬ 
ever,  recent  experimental  observations  [17,18]  sig¬ 


Elementary  Block  p* 


Fig.  4.  Physical  representation  of  shear  crack  analogue. 


nify  that  first  grain  boundaries  do  not  always 
cause  drops  of  growth  rate  as  a  small  fatigue 
crack  penetrates  it.  Some  grain  boundaries  take 
less  effect  of  blockage  to  small  fatigue  crack 
growth.  In  other  words,  Eq.  (12)  may  be  valid  for 
small  fatigue  crack  propagation  in  the  first  few 
grains.  The  Stage  I  growth  is  considered  to  termi¬ 
nate  at  a  transition  crack  length  2a0  at  which  the 
shear  mechanism  may  shift  to  tension  mecha¬ 
nism.  Apparently,  the  prediction  of  growth  rate 
in  Stage  I  entails  a  determination  of  shear  strain 
at  the  crack  tip  due  to  the  gliding.  Despite  of 
different  damage  mechanisms  possibly  associated 
with  the  Stage  I  growth,  we  still  adopt  the  con¬ 
cept  of  the  process  zone  to  configure  the  small 
fatigue  crack  tip  as  shown  in  Fig.  4.  A  growing 
small  fatigue  crack  in  this  case  is  thus  physically 
modelled  as  a  micronotch  with  a  tip  radius  of 
( d*)Q  that  is  supposed  to  be  equivalent  to  the 
process  zone  size  associated  with  the  transition 
crack  length  2ac  and  determined  using  Eq.  (8). 
Applying  Eq.  (7)  for  the  strain  distribution  in 
front  of  a  micronotch,  the  cyclic  shear  strain  can 
be  approximately  estimated  in  the  form  of 


A8(d*) 
2p*w(x- a)  ’ 


(14) 


where  p*  is  still  specified  as  the  microstruc- 
turally-affected-zone  size.  Noting  At  =  jAa 
sin[2(-n-/2)  -  (p)]  (<p  =  45°),  the  cyclic  shear  strain 
at  x  =  a  +  (d  *  )„  can  be  rewritten  in  the  form  of 

ir  <1  --9  »n[2(|  -»)] 

t - 2 t^e - 

(15) 
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A  maximum  shear  strain  theory  [19]  is  used  to 
estimate  strain-controlled  fatigue  in  the  process 
zone  and  measure  the  fatigue  resistance  to  maxi¬ 
mum  shear  strain, 

=  1.30^(2  N{)b  +  1.50e'((2  Nf)c.  (16) 

The  valid  use  of  Eq.  (16)  relies  on  the  prerequi¬ 
site  that  shear  strain  makes  a  major  contribution 
to  fatigue  damage.  By  Eq.  (16)  the  small  crack 
growth  rate  is  therefore  simulated  in  terms  of 
fatigue  damage  caused  by  the  shear  stress  in  the 
process  zone.  As  a  result,  the  growth  rate  law 
equation  for  a  single  small  fatigue  crack  in  shear 
mode  takes  the  form  of 

da  It  K' 

dN  -  2(^*  )o|!  1.95(o-f'  -  (Tm)e'f ) 


X 


2 


—  (1  +n')/b+c 


(17) 


Eqs.  (11)  and  (17)  provide  a  unified  model  to 
predict  growth  rate  of  small  fatigue  cracks  for 
both  Stage  II  and  Stage  I  propagation. 


5.  Predictions  and  discussions 


influence  of  local  microstructure.  Recent  study 
has  shown  that  kinetic  transition  from  Stage  I  to 
Stage  II  growth  along  with  the  microstructurally- 
affected-zone  size  can  be  analytically  simulated 
[20],  Improved  predictions  of  growth  rate  da/d N 
by  the  model  are  illustrated  in  Fig.  6  for  7075-T6 
and  7150-T6  aluminium  alloys. 

The  present  model  can  predict  the  growth  rate 
of  physically  small  cracks  completely,  using 
macroscopic  low-cycle  fatigue  properties  via  the 
process  zone  in  conjunction  with  local  microstruc¬ 
ture,  without  any  use  of  regression  of  experimen¬ 
tal  data.  Among  those  fatigue  properties 
in',  K',  Of',  e'f,  b,  c )  of  aluminium  alloys,  the 
cyclic  strain  hardening  exponent  ri  is  typical  be¬ 
cause  the  change  of  growth  rate  is  mainly  depen¬ 
dent  upon  it.  As  the  ri  varies  from  unity  to  zero  a 
polycrystalline  material  exhibits  from  a  perfect 
plasticity  to  a  perfect  elasticity.  However,  it  is 
hard  to  discriminate  each  individual  effect  of 
those  physical  quantities  on  small  crack  growth 
rate. 

As  described  earlier,  local  microstructural  di¬ 
mension/elementary  block  has  a  significant  influ¬ 
ence  over  a  growing  small  crack.  Grain  size  is 
usually  considered  a  factor  to  influence  small 
crack  growth.  As  a  first  step  in  modelling  an 
average  grain  size  is  used  to  measure  the  mi- 


in  actual  predictions  of  growth  rates  average 
grain  size  is  initially  used  to  represent  the  mi- 
crostructurally-affected-zone  size  p*  to  roughly 
take  the  microstructure  effect  into  account.  Re¬ 
gardless  of  applied  stress  levels,  the  stage  I  growth 
is  considered  to  terminate  at  the  first  grain 
boundary  and  the  transition  from  shear  mecha¬ 
nism  to  tension  mechanism  is  thought  to  start 
there.  Fig.  5  illustrate  a  typical  prediction  of 
growth  rate  by  using  Eqs.  (11)  and  (17),  which  are 
generally  reasonable.  However,  in  the  transition 
regime  at  lower  stress  level,  Eq.  (17)  predicts  a 
constant  transition  crack  length  at  which  the 
growth  rate  is  a  minimum  without  stress  effect. 
This  is  in  contrast  to  physical  intuition.  Obvi¬ 
ously,  instead  of  average  grain  size,  the  mi- 
crostructurally-affected-zone  size  p  *  is  a  key  pa¬ 
rameter  to  be  kinetically  determined.  It  is  possi¬ 
ble  to  take  a  proper  measure  of  the  extent  of  the 


Surface  Crack  Length,  2a  (mm) 

Fig.  5.  Predictions  of  growth  rates  of  small  fatigue  cracks  by 
the  developed  model. 
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crostructure  effect.  Hence,  the  growth  rate  is  crack  is  primarily  controlled  by  the  internal 

directly  associated  with  grain  size.  Experimental  strength  of  grain  boundaries  together  with  local 

investigations  to  7091  and  7050  [20]  aluminium  metallurgical  features  around  the  crack  tip.  Ei- 

alloys  have  revealed  that  large  grain  did  not  re-  ther  grain  boundary  blocking  or  slip  band  orien- 

sult  in  fast  growth  of  small  cracks.  A  comparison  tation  may  cause  decelerated  (or  accelerated) 

of  small  crack  growth  rates  for  2024,  2090,  6080,  growth  behaviour.  Other  obstacles  bounded  in 

7075,  7150,  7017  aluminium  alloys  with  their  grain  matrix,  such  as  inclusions  or  grain  boundaries  at 

sizes  did  not  show  a  clear  evidence  of  grain  size  subsurface,  may  also  interfere  with  small  crack 

dependence  of  small  fatigue  crack  growth  rate.  It  growth  causing  deflection  of  the  crack  path  or 

appears  that  the  growth  rate  of  small  fatigue  reducing  growth  rate.  Therefore,  it  seems  inade- 


Surface  Crack  Length  2a  (mm)  Surface  Crack  Length  2a  (mm) 

Fig.  6.  Comparisons  of  predicted  growth  rates  with  experimental  results  [15]  (a)  7075-T6A,  (b)  7075-T6B,  (c)  7075-T6C  and  (d) 
7150-T6E. 
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quate  to  correlate  the  microstructural  dimension 
p  *  to  an  average  grain  size.  Moreover,  examina¬ 
tion  of  measured  data  showed  that  the  decrease 
of  growth  rate  at  first  grain  boundary  was  not 
always  observed  in  experiment  as  it  was  expected 
by  the  grain  boundary  blockage  mechanism.  In¬ 
stead,  any  grain  boundary  besides  the  first  grain 
boundary  may  act  as  an  effective  barrier  to 
severely  deter  propagation  of  small  fatigue  crack 
before  the  Stage  I  growth  terminates.  In  this  case 
a  pile-up  behaviour  of  dislocations  ahead  of  a 
small  fatigue  crack  tip  as  a  result  of  blocking 
effect  of  grain  boundary  prevails  if  a  small  fatigue 
crack  stems  in  a  grain  boundary.  The  blocking 
effect  is  particularly  evident  in  the  transition 
regime  from  Stage  I  to  Stage  II.  The  microstruc- 
turally-affected-zone  size  p*  should  be  deter¬ 
mined  by  the  transition  condition  rather  than 
takes  the  average  grain  size.  On  the  other  hand, 
stress  level  is  well  known  to  influence  small  fa¬ 
tigue  crack  growth,  particularly  the  transition 
from  shear  mechanism  to  tension  mechanism. 

t  t  t  t 


y  t 


Fig.  7.  (a)  An  elliptical  crack  in  an  infinite  plate  subjected  to 
tension  loading;  (b)  conformal  mapping  of  an  ellipse  to  a  unit 


Shear  mechanism  acting  as  a  dominant  mecha¬ 
nism  in  Stage  I  propagation  takes  long  to  con¬ 
clude  at  lower  stress  level.  This  fact  has  been 
taken  in  accounted  by  means  of  defined  process 
zone  because  the  process  zone  size  is  strongly 
dependent  upon  stress  level,  plastic  zone  size  and 
microstructurally-affected-zone  size. 


6.  Conclusions 

Early  stage  of  small  fatigue  crack  growth  ap¬ 
pears  to  be  controlled  to  a  great  extent  by  mi- 
crostructurally-affected-zone.  The  microstruc- 
turally-affected-zone  acts  as  a  bridge  in  the  mod¬ 
elling  across  a  gap  between  micro-fatigue  be¬ 
haviour  and  macro-fatigue  behaviour  although  it 
is  initially  defined  by  taking  an  average  measure 
of  local  microstructural  features.  Secondly,  crack 
growth  rate  is  predicted  to  increase  with  increase 
of  stress  level  because  high  stress  level  causes  a 
large  plastic  zone  ahead  of  a  small  crack  tip.  As  a 
result  shear  crack  growth  may  terminate  at  a 
smaller  crack  length  if  applied  stress  level  is 
higher.  Contrarily,  the  period  of  shear  crack 
growth  may  last  long  at  lower  stress  levels.  How¬ 
ever,  since  the  process  zone  size  is  determined 
either  by  plastic  deformation  or  by  local  mi¬ 
crostructure,  it  is  still  hard  to  discriminate  which 
one  of  the  two  takes  priority  in  early  stage  of 
propagation. 

Using  microstructurally-affected-zone  and  pro¬ 
cess  zone  concepts,  a  small  fatigue  crack  tip  is 
logically  represented.  Based  on  experimental  in¬ 
vestigations  the  evolution  of  plastic  zone  associ¬ 
ated  with  a  small  fatigue  crack  is  quantitatively 
evaluated  which  follows  a  decreasing  pattern  of 
power  function  during  the  small  fatigue  crack 
growth.  The  evolution  of  the  normalised  process 
zone  size  d* /p*  is  thus  simulated  in  terms  of  the 
normalised  plastic  zone  size  rp/a  via  the  stress 
level  in  conjunction  with  the  local  microstruc- 
turally-affected-zone.  As  a  result  of  modelling  a 
quantitatively  unified  methodology  for  predicting 
growth  rates  of  small  fatigue  cracks  has  been 
developed  whose  pronounced  feature  is  that  it 
reveals  a  coherent  relation  connecting  small  fa¬ 
tigue  crack  growth  with  bulk  fatigue  properties  of 
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aluminium  alloys,  in  other  words,  connecting  the 
microscopic  fatigue  growth  to  the  macroscopic 
bulk  fatigue  behaviour  via  the  microstructurally- 
affected-zone.  However,  it  seems  inadequate  to 
refer  the  microstructurally-affected-zone  size  p* 
to  any  specific  metallurgical  size  of  aluminium 
alloys.  The  microstructurally-affected-zone  size  is 
related  to  a  microstructurally-related  condition 
that  determines  a  transition  from  Stage  I  growth 
to  Stage  II  growth. 

Appendix  A 

Crack  opening  displacement  derivation 

Consider  a  quasi-ellipse  contained  in  an  infi¬ 
nite  elastic  body  subjected  to  uniform  uniaxial 
tensile  stress  (see  Fig.  7(a)).  The  displacement 
function  for  this  infinite  body  with  an  elliptical 
crack  under  uniaxial  tension  a  is  given  by  the 
relation  [3] 

2p(u  +  iu)  =  k<£(z)  -<f/{z)-tft(z),  (18) 

where  u  is  the  displacement  in  x-direction;  v  is 
the  displacement  in  y-direction;  p.  =  E/ 2(1  +  v); 
k  =  (3  —  v)/(l  +  v)  in  plane  stress  and  k  =  3  —  4v 
in  plane  strain;  <£(z)  and  i friz)  are  the  analytical 
complex  functions.  The  solution  of  Eq.  (18)  is 
obtained  by  method  of  conformal  mapping,  trans¬ 
ferring  the  ellipse  in  z-plane  into  a  circle  of  unit 
radius  in  f-plane  (see  Fig.  7(b))  by  the  transfor- 


respectively.  All  complex  functions,  HO,  4>  (£), 
i !>({;)  and  <f>'(&/6>'(0,  can  be  determined  by  tak¬ 
ing  |  =  e'9  and  substituting  them  into  Eq.  (19). 
The  values  of  displacements  can  be  calculated  by 
separating  the  real  part  from  imaginary  part  in 
Eq.  (19)  that  is  substituted.  The  imaginary  parts 
in  Eq.  (19)  is  the  displacement  in  y-direction. 
This  displacement  is  expressed  as, 


(20) 


where  6  is  the  angle  between  any  point  of  inter¬ 
est  along  the  ellipse  face  and  the  major  axis  of 
the  ellipse.  The  crack  opening  displacement,  8,  at 
any  point  of  interest  takes  the  value  of  2v.  Under 
plane  strain  condition  (k  =  3  -  4v),  thus  8  can  be 
analytically  expressed  in  the  form  of 


(21) 


Since  crack  is  small  it  is  reasonable  to  assume 
that  this  ellipse  can  be  regarded  as  a  “line  crack”. 
Thus  c  becomes  the  semi-crack  length,  a  and 
b-*0,  consequently  resulting  in  M  =  1  and  A  = 
a/2.  After  these  simplifications  plus  x  =  a  cos  0, 
sin  0  =  [1  -  (x/a)2]l/2  in  (x,  y)  coordinate  sys¬ 
tem,  the  crack  opening  displacement  of  this  el¬ 
lipse  under  plane  strain  condition  (k  =  3  -  Av) 
can  be  calculated  by 


mation 

z-»(f)-^(f+ j)  M<  1, 

where  A  =  (c  +  b)/2  and  M  =  (c  -  b)/(c  +  b);  c 
and  b  being  the  semi-major  and  semi-minor  axis 
of  the  ellipse  respectively.  After  necessary  trans¬ 
formation,  Eq.  (18)  becomes 
2 n(u  + iv) 


Eq.  (22)  is  similar  with  Westergaard’s  crack  open¬ 
ing  displacement  equation.  When  plastic  zone  of 
size  rp,  that  is  a  distance  from  the  crack  tip  to  the 
edge  of  the  plastic  zone,  is  formed  at  this  ellipti¬ 
cal  crack  tip  we  considered  an  equivalent  elastic 
elliptical  crack  whose  tips  in  x-direction  are  lo- 


- k<ho  -  nowm^m  -  m-  m 

Above  equality  can  be  solved  by  letting  [22] 


HO 


2  +  M\ 

~Tr 


HO 


<rA  U3-(M2  +  2M+1)€-Cl 
M  e-M 


cated  at  the  edge  of  the  plastic  zone.  Therefore, 
8  is  defined  as  the  distance  of  separation  of  two 
faces  of  this  equivalent  crack  due  to  its  notional 
growth.  The  corresponding  crack  opening  dis¬ 
placement  of  this  equivalent  crack  is  thus  calcu¬ 
lated  by  replacing  a  with  a  +  rp  in  Eq.  (22). 
However,  since  the  crack  is  very  small,  the  crack 
is  actually  embedded  in  a  local  plastically  de- 
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formed  region  where  the  hypothesis  of  small  scale 
yield  is  no  long  applicable.  Plenty  of  experimental 
research  on  notch-like  crack  [23-26]  suggests  that 
a  crack  mouth  opening  displacement  (CMOP) 
may  be  adequate  for  a  fully  elastic-plastic  condi¬ 
tion  to  represent  the  actual  crack  tip  opening 
displacement.  For  a  surface  small  fatigue  crack, 
take  x  -  a  in  Eq.  (22)  to  determine  the  crack  tip 
opening  displacement. 
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